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ABSTRACT 

Context. Recent measurement of mass of PSR J1614-2230 rules out most of existing models of equation of state (EOS) 
of dense matter with high-density softening due to hyperonization or a phase transition to quark matter or a boson 
condensate. 

Aims. We look for a solution of an apparent contradiction between the consequences stemming from up-to-date hyper- 
nuclear data, indicating appearance of hyperons at ~ 3 nuclear densities and existence of a M = 2.0 Mg neutron star. 
Methods. We consider a non-linear relativistic mean field (RMF) model involving baryon octet coupled to meson fields. 
An effective lagrangian includes quartic terms in the meson fields. The values of the parameters of the model are 
obtained by fitting semi-empirical parameters of nuclear matter at the saturation point, as well as potential wells for 
hyperons in nuclear matter and the strength of the A — A attraction in double-A hypernuclei. 

Results. We propose a non-linear RMF model which is consistent with up-to-date semiempirical nuclear and hyper- 
nuclear data and allows for neutron stars with hyperon cores and M > 2 M©. The model involves hidden-strangenes 
scalar and vector mesons, coupled to hyperons only, and quartic terms involving vector meson fields. 
Conclusions. Our EOS involving hyperons is stiffer than the corresponding nucleonic EOS (with hyperons artificially sup- 
pressed) above five nuclear densities. Required stiffening is generated by the quartic terms involving hidden-strangeness 
vector meson. 
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1. Introduction 

Recent measurement of the mass of PSR J1614-2230, 1.97± 
0.04 M Q (Demorest et al. 2010), puts stringent constraint 
on the equation of state (EOS) of dense matter in neutron 
star cores. In the light of this measurement, EOSs of dense 
matter, based on the modern many-body theories and real- 
istic strong-interaction model, lead to a puzzle. On the one 
hand, interactions consistent with experimental data on hy- 
pernuclei, lead to the presence of hyperons at the densities 
exceeding 2 — 3po, where po = 2.7 x 10 14 g cm -3 (corre- 
sponding to the baryon number density rig — 0.16 fm ) 
is the normal nuclear density. On the other hand, in- 
evitable softening of the EOS, due to the hyperonization, 
imp lies the maximum allowable mass M max ^1-5 Mq (see, 
e.g. lBurgio et al. ll201lllVidana et al. 11201 lL and references 
therein). Such a low M max is only marginally consistent 
with 1.44 M Q of the Hulsc- Taylor pulsar, but was con- 
tradicted already by 1.67 ± 0.04 M of PSR J1903-0327 
(Champion et al. 2007; more precise value has been re- 
cently obtained by Freire et al. 2010). Let us mention, that 
this problem cannot be solved by adding an ad hoc ex- 
tremely stiff repulsive three-body contribution to the EOS 
(jVidana et al. 1 . 120 111 ). 

Wc consider neutron star cores composed of baryons, 
electrons, and muons. Baryons and leptons are in weak- 
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interaction equilibrium. Hyperons appear at density p± 
(baryon density n{). For p < pi only nucleons are present. 
For p > p\ matter contains a mixture of nucleons and hy- 
perons. This state (phase) will be denoted as NH. We will 
also consider a dense matter model with hyperons artifi- 
cially suppressed. Such a purely nuclcon state will be de- 
noted as N. The corresponding EOS will be denoted as 
EOS.NH and EOS.N. These EOS coincide for p < p x . 

A too low Min^ is not an inevitable feature of neu- 
tron stars with hyperon cores. In a very recent paper 

iBonanno fc Sedrakian I (|2011l) get M§S> > 2 M Q , start- 
ing from an extremely stiff relativistic mean field model 

NL3 EOS.N, which yields M^ x = 2.8 M© (close to the ab- 
solute upper bound o n M max stemming from causal ity, e.g., 
lHaensel et"aT1l2007t) . fBonanno fc Sedrakian I (|2011l) extend 

the NL3 model to the hyperon sector, and get AftSx = 
2.03 Mq. Massive stars with hyperon core exist there be- 
cause of the extreme stiffness of EOS.N, obtained at the 
cost of having a parameter of nuclear-matter parameter 
(the slope of symmetry energy vs. density) L = 118 MeV 
significantly higher than its semi-empirical estimates. 

In our approach we keep L (and the stiffness of the 
EOS.N near po) within the semi-empirical (i.e., obtained 
within a model of atomic nuclei, and therefore model- 
dependent) estimates. Our EOS.N at high density is stiff, 

but not extremely stiff: M^X = 2.1 M . Generally, Mmax' 
is essentially determined by the p 5po segment of the 
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EOS.NH. So if the hyperon softening occurs at 2 — 3po, 
then, to get a sufficiently large Af max , the softening must 
be followed by a sufficiently strong stiffening of EOS.NH for 
p <; 5po- One has therefore to find a mechanism of stiffening 
of EOS.NH at these densities. For our model of EOS.N, in 

order to yield m£S? > 2 M Q , EOS.NH for p £ 5p should 
be actually stiffer than the EOS.N one. Simultaneously, the 
NH phase has to be stable (thermodynamically preferred 
over the N one). We derive a constraint on the EOS.NH 
resulting from the conditions mentioned above, and discuss 
consequences of violation of this constraint. 

In our discussion we restrict to hadronic matter, and do 
not consider a possibility of quark deconfincment. Because 
of the surface effects and electrical screening in a quark 
plasma, transition to a quark matter would occur at nearly 
constant pressure and with an only slig htly smoothed 
density jump (see, e.g.. lEndo et al. 1120061 and references 
therein). A reasonable scenario is: softening of the N phase 
by hyperonization at 2 — 3po, then softening of the NH 
phase by quark deconfinement at a significantly higher 
density. iBonanno fc Sedrakian I ([20111 ) show that assuming 
NL3 EOS.NH, a transition to quark matter occurring after 
hyperonization could be consistent with 2 M Q pulsar pro- 
vided the vector repulsion in quark matter is sufficiently 
strong and quark deconfincment takes place near the maxi- 
mum NS mass. For our EOS.NH getting 2 M with a quark 
core would require a very fine tuning, and we do not con- 
sider such an unlikely possibility. 

The problem of an interplay of attraction (softening) 
and repulsion (stiffening) in dense hadronic matter can 
be formulated in simple terms using a modern effective 
field theory, involving baryon and meson fields. In the 
case of nucleon matter, such a theory can be p ut on the 
firm theoretical basis, starting from the QCD ([Walecka I 
12004 and references therein). Such an effective theory 
can be solvable within the mean field approximation and 
can give a satisfactory description of a wealth of nuclear 
physics data if the coupling of nucleons to the three me- 
son fields: scalar a, vector w M and p 1 ^ , is considered. Here, 
p and i denote the space-time and isospin-space compo- 
nents of the field. Effective lagrangian contains quadratic 
and quartic terms in vector fields, and quadratic, cubic, 
and quartic terms in scalar fields. While a yields attrac- 
tion to bind nuclei, vector meson fields generate repul- 
sion to saturate nuclear matter at po. Numerical coeffi- 
cients in the effective l agrangian are fixed by fitting a 
wealth of nuclear data (|Sugahara fc Toki I Il994f) . The ef- 
fective model is then extended to include the hyperon sec- 
tor. Two meson fields with "hidden strangeness" (ss) are 
added: scalar a* (quadratic terms) and vector-isovector <j) % 
(quadrati c and quartic terms). T hese fields couple to hyper- 
ons only (Schaffn er et al. II1994I ). An important constraint 
on the hyperon sector of lagrangian results from the ex- 
isting evaluations of the depth of the potential well act- 
ing on a single zero-momentum hyperon in nuclear mat- 
ter, u[ N \ C/j, ,and f/n (binding energy of a hyperon 



in nuclear matter is Bh 



— U H ). An effective theory 
of hadronic matter described in general terms above and 
solved in the mean field approximation will be referred to as 
non-linear relativistic me an field model (non- linear RMF, 
iBednarek fc Mahkal [20091 ). 

There exist a few older and simpler models of NH 



aim by pushing the parameters of the interaction to the 
extreme and likely unrealistic limits. These models are 
based on the relativistic mean field lagrangian involving 
octet of baryons coupled to a, lu^ and p 1 ^ meson fields 
(quadratic terms in the lagrangian), with additional cu- 
bic and quartic a s elf-interaction terms (for a review, see 
IGlendenning i n"99& i. The mean field solutions of the field 
equations are called relativis t ic me an field (RMF) model. 
iGlendenning fc Moszkowski I (|l99lh can exceed 2 M Q as- 
suming unrealistically high nuclear matter incompressibil- 
ity, K = 300 MeV, and an unrealistically strong A — a 
attraction, balanced by a A — w repulsion to get exper- 
imental Uj^\ Relations between coupling constants gnu> 
and gNta, resulting from the SU(6) symmetry and making 
them significantly weaker than g^u, have to be strongly 
violated. Finally, experimental constraints on 17s and Us 
in nuclear matter are not applied. Instead, it is assumed 
that all hyperons in the baryon octet have the same cou- 



(NH) 
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pling as A. Similar mode l was used to obtain M, 
2 M bv iBombaci et al. I ([2008D . Very recently, the RMF 

model was applied to calculate A/„iax ' under the conditions 



U 



<iV> 



-30 MeV uj N) = +30 MeV, U±" > = -18 MeV 
and varying gnuj (jChatteriee fc Schaeffner-Bielich I [201 it ). 
For the SU(6)-symmetric set of gnu, the maximum mass 
is about 1.7 M Q . Mmlx > 2 M Q can be obtained only for 
sufficiently strong giiu — SNlo, i-e. by breaking the SU(6) 
symmetry. 

In the present paper we propose a solution of the puz- 
zle "hyperonization — M max > 2 M Q " using a specific re- 
alization of the non-linear RMF model of hadronic matter 
([Bednarek fc Mahkal I2009T) . The non-linear RMF model of 
NH matter is presented in Sectf2] Experimental constraints 
from nuclear and hypernuclear physics are described in 
Sect 01 A particular EOS.NH is described in SectH A 
model of PSR J1614-2230 with a hyperon core is presented 
in SectJSJ The problem of high-density instability of the NH 
phase and the M — R relation for neutron stars models are 
discussed in SectEl Finally, in Sectf7]we summarize results 
of the paper, compare them with results obtained by other 
authors, and present our conclusions. 

Preliminary results of our work were presented at 
the MODE-SNR-PWN Workshop in Bordeaux, France, 
November 15-17, 2010, and in a poster at the CompStar 
2011 Workshop in Catania, Italy, May 9-12, 2011. 



2. Non-linear RMF model of hyperon cores 



matter, yielding M, 



(NH) 



> 2 Mq, but they reach this 



The m odel has been formulated by IBednarek fc Mahkal 
([2009( 1. The octet of baryons includes nucleon doublet N 
and six lowest-mass hyperons H: A singlet, S triplet and S 
doublet. Uniform number density of each baryon species B 
is denoted ub (B = n,p, A, . . .). 

In the nucleon sector the meson fields arc: scalar a, vec- 
tor uj^, and vector-isovector p 1 ^. Generalization of the non- 
linear RMF model to the baryon octet is done in the follow- 
ing way. Additional "hidden-strangeness" mesons: scalar a* 
and vector c/> M are introduced. They couple to hyperons 
only, gNa* = 9Nc/> = 0. The vector-meson coupling con- 
stants to hyperons are assumed to fulfill relations stemming 
from the SU(6) symmetry (additive quark model): 
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Table 1. BM165 model of hadronic matter. Calculated nu- 
clear matter parameters at the saturation point 



3Aw = .9Sw = 2.g E 



^gNw , gup = 2gs P = 2g Np , (l) 

gs</> V2 
gA<t> — .9s0 — —j- — —gNu- 





E s 


K B 




L 


(fin" 3 ) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


0.145 


-16.3 


279 


33 


74 



Similar symmetry relations can be obtained for the coupling 
constants of the scalar mesons, but they are not used in 
the present model. Instead, we are adjusting them to fit 

experimental estimates of U B N \ 

At fixed {n-s} (B = n,p, A, . . . ), and assuming vanish- 
ing baryon currents, the hadronic lagrangian density Chad 
is used to derive equations of motion for the meson fields. 
Static solutions are found assuming that baryonic matter is 
isotropic and uniform. Mean-field approximation, neglect- 
ing quantum corrections, is used: o — > (a) = sq, lu p — > 
(ujp) = w Q 5 pQ , — > (p'l) = ro^ofe, o* — > (a*) = s* , 
<f> p — > (<f> p ) = foSpQ. The resulting lagrangian density 
function £had consists of three components, £had = £b + 



C^' +C]^'' 1 ' . Cb is obtained from the free-baryon lagrangian 

,...) 
The 



,(3,4) 



by replacing bare baryon masses mj (B = n,p,A 
by the effective ones, m*. 



quadratic (interaction) component £^ contains terms pro 



tub - gBaSa - gBa* s 

(2) 



portional to Sq, Wq, Tq, Sq , /q . The interaction component 
£^' 4 ^ contains terms cubic and quartic in sq, and quartic 
vector- meson terms proportional to Wq, rfi, /q and the cross 
terms proportional to fnWn, f 2 rf, and 



wZr, 



o'o- 



The hadronic lagrangian density function £had is then 
used to calculate the hadron energy-density as a function 
of partial baryon densities (tibI , 4a.H (I^bD- Calculations 
done for the considered model (|Bednarek fe Mahka 20091 ) 
lead to the following explicit formulae for the hadron con- 
tribution to the energy d ensity £ and pressure P (n otice 
that original equations in iBednarek fc Mahkal (|2009l) con- 
tain several misprints which are corrected below; for the 
sake of simplicity, we use a shorthand notation g^ a = 

g<j , gNu = guj , • • ■ ) 



^had = +\mlsl + \mlwl + \m 2 p rl + U(s ) 



y- 

^2 

B 



2 

kF,B 



k dkJ fc 2 + (m B - g„BSo 



-gv s o 



+3A v (g p g u ) 2 w 2 r 2 + -c 3 (wf, + r^) 



,1 2 f 2 , 1 2 *2 



I lc 3 + ^A v (g p9u ,y 



( \c 3 - l -K v (g pgbJ f 



(2) 



fhad = -^mlsl + ^m p r 2 + ^m^w 2 -U(s ) 



y- 

^ 3tt 2 



fc 4 dfc 



(fc 2 + m B - gaBSo 



-g<?- s o 



+-c 3 (^ + t-q) + Kv{g P gLo) 2 wlrl 



i 1 2 j-2 1 2 *2 



-c 3 +-A v (g pgui y 



[ \^v{g P g^f ) JoVo + '-o) 



(3) 



where the non-linear cr-sclf-intcraction potential is 

U(°) = -^-T^- (4) 



The terms vanishing in purely nucleon (zero strangeness) 
matter are put into the rectangles. 

The quartic terms in f had and Phad deserve additional 
explanations. Their form stems from the chiral SU(3) sym- 
metry of the baryon-meson and meson-meson interactions. 
The coefficients of the quartic terms involve two phe- 
nomcnological parameters, c 3 and Ay. 

Let us consider first the quartic terms in the nucleon 
sector. The vector-isoscalar qua rtic term (uyj) was alread y 
present in the TM1 model of ISugahara fc Toki I (|1994[ ). 
However, TM1 was constructed to describe atomic nuclei, 
and therefore it is valid for nuclear matter near saturation 
density and for small neutron excess. It is to be expected 
that extrapolation to supranuclear density and large neu- 
tron excess necessitates a richer isospin and density de- 
pendence of the model lagrangian than t hat assumed in 
the TM1 model. [Bednarek fc Mahkal (|2009h proposed to do 
this by enlarging the quartic terms via adding a vector- 
isovector one (r-g) and a cross-term (w 2 r 2 ). The strengths 
of the quartic terms are determined by two parameters, 
Ay and c 3 , instead of only one in the TM1 model of 
ISugahara fc Tokil (|1994|) . This allows for a good fitting not 
only (semi-empirical estimates of) nuclear symmetry en- 
ergy and incompressibility, but also the slope parameter L, 
and simultaneous ly yields M max > 1.97 M . 

As shown in (jBednarek fc Mahkal l200l . chiral SU(3) 
symmetry yields a suitable extension of the quartic terms 
to the hyperon sector, the same c 3 and Ay entering the 
quartic-terms coefficients. Additional quartic terms in the 
hyperon sector are generated by the hidden-strangeness 
vector-isoscalar field /q. 
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Table 2. BM165 model of hadronic matter. Calculated 
zero-momentum single-baryon potentials in symmetric nu- 



clear matter, 



and of A in A matter, U A 
are obtained at nb = n s 



(A) 



All results 



U (N) 

N 

(MeV) 


U A 

(MeV) 


(MeV) 


(MeV) 


U A 

(MeV) 


-69 


-28 


+30 


-18 


-5 



3. Determination of parameters of non-linear RMF 
model 

Let us denote the neutron excess in nuclear matter by 
8 = (n n — rip)/ rib . The energy per nucleon (excluding nu- 
clcon rest energy) is E(nb, 8). Analysis of a wealth of data 
on heavy atomic nuclei can yield the parameters of nu- 
clear matter near the saturation point, corresponding to 
the minimum of energy per nucleon, E s , reached at = n s 
and 8 = 0. Results are mo del- dependent and therefore they 
are called semi-empirical. Other semi-empirical parameters 
arc: symmetry energy S s , incompressibility K s , and the 
symmetry-energy slope parameter L, 



S a = 



d 2 E 
d8 2 



= 3n = I - 



n a ,8=0 

d 3 E \ 



f d 2 E 



,5=0 



\9n h d8 2 J 4=0 



(5) 



Studies of hypernuclei and of E atoms allow for eval- 
uation of the potential energy of a single zero momentum 

hyperon in symmetric nuclear matter, f/jy . The non-linear 
RMF yields following expression for this quantity: 



U 



(N) 
H 



9HcrS0 - 9HujW , 



(6) 



4 S) and l4 H) 
Therefore, we estimate them using the relations 



used to determine g\&* ■ For and no data exist. 



U, 



(2) 



u 



(2) 



2U. 



(A) 



2U 



(A) 



(8) 



These relations have been established on the basis of one- 
boson exchange models and semi- empirical evaluation o f 
the strength of the A — A attraction (jSchaffner et al. II1994I ). 

We adjusted parameters of our lagrangian to repro- 
duce, within a few percent, ten semi-empirical nuclear and 
hyper-nuclear data. Adjusting the isovector parameters in 
the lagrangian density, g p and Ay, deserves an additional 
explanation. This was done by employing existing infor- 
mation on the density dependence of the symmetry en- 
ergy. Namely, we used not only the value of the symme- 
try energy at saturation, S s , but also a semi-empirical es- 
timate of symmetry energy at rib ~ 0-1 fm~ 3 , 26.67 MeV, 
which determined the value of L (|Horowitz fc Piekarewiczl 
l200ll ). This influenced the EOS of neutro n matter, be- 
cause E{rib, 1) ~ E{rib, 0.5)+5(nb) (see, e.g.. lHaensel et all 
[2001 
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Fig. 1. Equations of state EOS.N and EOS.NH calculated 
using the BM165 model. Hyperons appear at point 1 and 
EOS.NH crosses the EOS.N one at point 2. 



which should be calculated at n s and 8 = 0. The scmi- 



(iV) 
A 



-28 MeV, U. 



empirical estimates are U 

+30 MeV, ui N) = -18 MeV (|Schaffner-Bielich fc Gal I 
120001 ). Eq. ((5J) is then used to determine g\_ a , g^ a and g-^a- 
As we are dealing with NH phase, containing finite frac- 
tions of hyperons, we need information on the hypcron- 
hyperon interaction. Studies of double-A hypernuclei sug- 
gests that A — A interaction is attractive. In the mean-field 
approximation it can be characterized by the potential well 
of a zero momentum A in the A-matter. In our model we 
get a general expression for a potential energy of a zero- 
momentum hyperon H' in the //-matter, 



in 1 \ 

U H = gHaSo - 9HluW + 9H<j* Sq - gH<t>fo ■ (7) 
The latest (very uncertain) semi-empirical estimate 



coming from double-A hypernuclei is U 



(A) 



-5 MeV 



(jTakahashi et al. II2001I . fSong et al. II2003D . Eq. © is then 



4. EOS of neutron-star matter 

Total energy density and total pressure are sums of con- 
tributions of hadrons and leptons, because contributions of 
electromagnetic interaction to these quantities are negligi- 
bly small, and leptons (electrons and muons) can be treated 
as ideal Fermi gases, 



£ — £had + £lcp , P — PhaA + Plep 

Let us fix the (total) baryon number density, 



nb 



E 



II B 



(9) 



(10) 



Imposing electrical charge neutrality and beta equilibrium, 
one gets a system of non-linear equations for the particle 
species fractions Yi = n^/rib- Let at a given nb the equi- 
librium fractions be {Y i (rib)}- Under imposed conditions, 
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Fig. 2. Number fractions of the constituents of dense mat- 
ter in beta equilibrium, Yj = rij/rib, versus baryon density, 
n b . Dotted lines: EOS.N. Solid lines: EOS.NH. 



energy density and pressure become functions of rib only. 
For the sake of comparison, wc will consider not only a gen- 
eral EOS involving nucleons and hyperons, EOS.NH, but 
also an EOS of nucleon matter, EOS.N, with artificially 
suppressed hyperons (this was done by doubling the actual 
hyperon rest masses). 

Expression for Phad (Eq[3]) shows that hyperons produce 
new repulsive quartic terms involving /o, u>o, and r$. The 
dependence of the EOS on Ay turns out to be quite strong. 
All other parameters being determined by the conditions of 
reproducing nuclear and hypcrnuclear data (i.e., near the 
saturation point of nuclear matter) , the high-density stiff- 
ness of EOS.NH increases monotonically with increasing 
Ay. In what follows we will use Ay = 0.0165, which as 
we will see yields a high-density stiffness of the EOS.NH 
consistent with M„Hx' > 2 M Q , while keeping good agree- 
ment with semi-empirical nuclear matter parameters. This 
EOS.NH will be hereafter referred to as BM165.NH or 
BM165. Thermodynamical equilibrium in the NH matter 
imposes relations between chemical potentials of hyper- 
ons present in dense matter, nucleons, and leptons. These 
relations can be ex pressed in the general form (see, e.g., 
lHaensel et aTll2007D : 



He = 



(11) 



where q H is the charge of hyperon H in units of the proton 
charge. The threshold density for appearance of hyperons 
H , , is determined by the enthalpy of a single hyperon 
H in a bcta-equilibratcd dense matter, 



/% = m H c 2 + U H + P/n h 



where Uh is the potential energy of this hyperon. For rib < 
> fi n —q H fjt & , and H decays in a weak interaction 
process. At the threshold density a single H in dense matter 
is stable, 

/% = Mn -1 H ^e , (13) 

and for rib > the density of stable H grows with in- 
creasing nb- 



For the BM165 model first hyperon to appear is A at 
6.3 x 10 14 g cm -3 (0.35 fm -3 ), where ^ = /i„. Next hy- 
peron is S", it appears at 7.7 x 10 14 g cm" 3 (0.42 fm -3 



), where = \i n + /i e . Repulsive U^l' implies that S 
appears at significantly higher density than S~, namely 
1.23 x 10 15 g cm -3 (0.63 fm -3 ), in spite of m s - < m 3 -. 

The order of appearance of hyperons in dense mat- 
ter deserves a comment. For a long time, in view of lack 
of an experimental information on E/j. and lli N \ they 

were assumed to be similar to . Consequently, £~ was 
found to be the first hyperon to appear, not the lightest 
hyperon A, because the (unfavorable) effect of ms- > tba 
was weaker than the (favorable) effect of the presence of 

/ip in the threshold condition p,^l — fi n + /i e (see, e.g., 
lHaensel et all l2007t ). However, large positive (repulsive) 

resulting from analyses of the £ - -atoms pushes 

well above ri£ . Consequently, £~ is the last, instead of 
being the first, to appear in neutron star core. 

Appearance of hyperons leads to significant softening of 
the EOS.NH compared to EOS.N (Fig. 2). In order to sup- 
port neutron stars with M > 2 M Q , the EOS.NH has neces- 
sarily to significantly stiffen at higher densities. The curve 
P NH (p) crosses the P N (p) one at p 2 = 1.76 x 10 14 g cm -3 
(n 2 = 0.85 fm -3 ), and at higher density EOS.NH is stiffcr 
than EOS.N. 

Actually, the difference P NH — P N is limited by the con- 
dition of stability of the NH phase against the re-conversion 
into the N phase. Assume that matter is in beta equilib- 
rium. At T = the small change of energy per baryon dP 
is related to the small change of baryon density drib by 



dE = P 



dni. 



(14) 



Therefore, condition of stability of the NH phase (against 
the conversion into the N one) 



p nh M <p 1n K) 



implies 



- b p*H K) -P»K) dn , <0; 



(15) 



(16) 



where n% is the density at the first hyperon threshold. For 
the BM165 model n x = 0.35 fm -3 . 



5. A model of PSR J1614-2230 

To get a complete EOS of the neutron-star interior, the 
BM165 EOS of the liquid core was supplemented with an 
EOS of the crust. We use d the EOS o f the inner crust of 
iDouchin fc Haensel (|2001l) , the model of lHaensel &: Pichon I 



(12) (|1994[ ) for the ou ter crust down to 1 8 g cm 3 , and the 



classical model of iBavm et al. I ()1971l ) for the outer layer 
with p < 10 s g cm -3 . A model of neutron star of grav- 
itational mass 1.97 M Q , rotating rigidly at 317 Hz, was 
calculated using the 2-D rotstar code from the LORENE 
library ( jhttp : //www . lorene . obspm . fr| implementing the 
formulation of (|Bonazzola et al . 1993). The circumferential 
equatorial radius of neutron star is P cq = 11.83 km, central 
density p c = 1.73 x 10 14 g cm -3 and central baryon den- 
sity n c = 0.834 fm -3 . At 317 Hz, polar flattening is rather 
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Fig. 3. EOS BM165 in the vicinity of the high-density NH- 
N phase transition. 2: crossing of the the N and NH pres- 
sures. 3: density and pressure of the NH phase at the phase 
coexistence. 3': density and pressure of the N phase at the 
phase coexistence. Horizontal segment 33' - pressure at the 
1-st order phase transition. Were the NH phase be (suf- 
ficiently) stable till the maximum mass, maximum (cen- 
tral) density in stable stars would correspond to the as- 
terisk sign, Pc,max = P*- For the fully equilibrated hadronic 
matter p Ciimx is slightly higher than p%. For a more detailed 
discussion of this point see the text. 




Fig. 4. The logarithm of the number fractions of the con- 
stituents of dense matter, log 10 (Yi), versus circumferential 
radius, in the liquid core of a 1.97 M© star model based on 
the BM165 EOS. 



small: the radial coordinate at the equator is only 200 m 
larger the that at the pole. The number fractions of the 
particle species Yi = n^/rib, versus radial coordinate r, in 
the liquid core of neutron star, are plotted in Fig. [4] The 
radius of the hyperon core is 8.36 km. The strangeness per 
baryon at the star's center (SjN^) = —0.35. 



6. High-density instability of the NH phase and 
neutron star models 



Violation of inequality (fT6"|) indicates instability of the NH 
phase against it conversion into the purely nucleon (N) one. 
Thermodynamic equilibrium of dense matter at pressure P 
corresponds to the minimum of baryon chemical potential 
/ib = {£ + P)/nh- An equilibrium phase-transition NH — >~N 
occurs at P3 such that 



/TVs) 



Mb 1^3 ) , 



(17) 



and is accompanied by a density jump from p% = 
p(NH)(-p 3 ) = 2.58 x 10 15 gem" 3 (n 3 = 1.105 fm" 3 ) on 



the NH side to p' z 



,(N) 



(Ps) 



3.25 x 10 i5 g cm 



_ (n' 3 = 1.31 fm" 3 ) on the N side. The BM165 EOS in the 

4xio 15 vicinity of P3 = 1.05 x 10 36 erg cm -3 is shown in Fig. [31 
The softening of the EOS for P > P 3 is twofold. First, 
there is a constant pressure sector of the EOS (vanishing 
compression modulus!). Second, there is a transition to the 
N-phasc which is significantly softer than the N one. 

The reaction of the star structure to the (1st order) 
phase transition (NH — s>N) can be described by the linear 
response theory formulated in IZdunik et al.l (|1987[ ) . This 
theory describes stellar configurations in the vicinity of the 
star with central pressure equal to P3. In our case this re- 
gion of stellar configurations is very small since we are close 
to the maximum mass. The crucial parameter, determining 
the stability of the star with a small core of the denser 
phase (N), is the relative density jump at phase transition 
pressure, P3: A = p' 3 / P3- Stability condition for a star with 
a small N-core reads: 



A < A crit = -(1 



P3C Z 



(18) 



(see Sect 3.4 of IZdunik et"al] [l987h . In our case condi- 
tion (TT5|) is fulfilled, because A = 1.27 while A cr it = 2.18. 
Consequently, there exists a (very small) region of stable 
configurations with the N-phase core. Actually, this region 
is very narrow: the maximum mass of non-rotating stars 
is reached for central pressure P c ,max larger only by 0.04% 
than the pressure at the phase transition, P3. 

The M — R relation for non- rotating NS, and those ro- 
tating at 317 Hz, based on the BM165 EOS, are plotted 
in Fig.[S] Stars with M > 1.4 M Q have a hyperon core. 
The flattening of the M{R) curve due to the hyperon soft- 
ening of the EOS is significant. However, it still allows for 
M ^NH)stat = 2 Q3 which ig Qnly by Q Q7 Mq smaller 

than M^lx tat . Rotation at 317 Hz, as measured for PSR 

J1614-2230, increases M, 
in Fig.0. 



(NH) 



to 2.04 Mq (see zoomed inset 



7. Discussion and conclusions 

We constructed a model of hyperon cores of neutron stars 
that allows for the existence of neutron star of 2 M . The 
model is consistent with ten semi-empirical evaluations of 
nuclear and hypcr-nuclear matter parameters. As an addi- 
tional constraint, we imposed SU(6) symmetry relations be- 
tween the coupling constants of baryons and vector mesons. 
In spite of this, by introducing two hidden-strangenes me- 
son fields (scalar and vector) coupled to hyperons only, 
we were able to reproduce four semi-empirical parameters 
stemming from hypcrnuclcar physics. 



I. Bednarek et aL: 2 Mq pulsar and hyperon cores 



7 




Fig. 5. Gravitational stellar mass, M, versus circumferen- 
tial radius, R, calculated for the EOS.N and EOS.NH. Only 
stable configurations are displayed . Inset: effect of rotation 
at / = 317 Hz on the M - equatorial circumferential radius 
curve near M max . 



In contrast t o NL 3 model, used by 
iBonanno fc Sedrakian I (poll , our symmetry energy 
is not unusually "stiff' near the saturation point: we 
;et L = 74 MeV, co mpared to L = 118 McV for NL3 
Agrawal et"al] l2005t) . Consistently, our EOS.N is not 
unusually stiff, and yields for NS with nucleon cores 

M,S2c 8tat = 2.10 Mq, to be contrasted with the NL3 
va lue of 2.8 M^. The hyperon softening for the model 



of IBonanno fe Sedrakian" (|201ll )) is dramatic, and leads 



to M 1 ^x' )stat which is lower by nearly 0.8 M Q than 
the N one. In our case, getting M^^ stat > 2.0 M© is 
conditioned by the high-density vector interactions in 
the hyperon sector (not excluded in view of our lack of 
knowledge of high-density hyperon interactions) , while 
IBonanno fc Sedrakian I (|201ll ) rely on the assumed extreme 
stiffness of the EOS.N, which is necessarily connected with 
a very high value of their L. 

Our EOS.NH becomes stiffer than EOS.N for p > 5/9 , 
and its stiffness grows with density. This leads eventually 
to the instability of the NH matter against the conver- 
sion into the N one, softening the EOS due to the first 
order phase transition. The maximum density at which sta- 
ble NH phase can exist determines actually our M^lx' stat , 



which is only 0.07 Mq lower than M^[lx tat . 



Rotation at 



317 Hz, measured for PSR J1614-2230, increases Afff by 
0.01 Mq, to 2.04 M . Breaking the SU(6) symmetry for 
the vector-me son couplings to hyperons, similar l y as i n the 
recent work of lChatteriee fc Schaeffner-Bielichl ((201 it ) , can 

make the value of MmJ>P significantly higher (Bednarek et 
al., in preparation). 

In the present paper we restricted to the hadronic 
forms of matter. A consistent treatment of the phase 
transition to quark phase in neutron star core would re- 
quire the use of the QCD for both hadronic and quark 
phases. As the transition occurs in the strong-coupling 
regime, approximations are not controllable. An approach 
based on an effective model of the QCD of quark mat- 



ter (Na mbu-Jona-Lasinio) and NL3 for hadronic phase, 
used bv IBonanno fc Sedrakian I (|201lh indicates that to get 
Mmax > 2 Mq vector repulsion in quark matter should 
be sufficiently strong. In any case, the maximum mass ob- 
tained by them is very close to that reached at central den- 
sity equal to the deconfinement density. 
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